Basins of attraction on random topography 
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We investigate the consequences of fluid flowing on a continuous surface upon the geometric and 
statistical distribution of the flow. We flnd that the ability of a surface to collect water by its mere 
geometrical shape is proportional to the curvature of the contour line divided by the local slope. 
Consequently, rivers tend to lie in locations of high curvature and flat slopes. Gaussian surfaces 
are introduced as a model of random topography. For Gaussian surfaces the relation between 
convergence and slope is obtained analytically. The convergence of flow lines correlates positively 
with drainage area, so that lower slopes are associated with larger basins. As a consequence, we 
explain the observed relation between the local slope of a landscape and the area of the drainage 
basin geometrically. To some extent, the slope-area relation comes about not because of fluvial 
erosion of the landscape, but because of the way rivers choose their path. Our results are supported 
by numerically generated surfaces as well as by real landscapes. 

PACS numbers: 92.40.Fb, 94.40.Gc, 05.90.+m 
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I. INTRODUCTION 

Aside from their natural beauty O the morphological 
features of landscapes also bear the signatures of tec- 
tonics and past and present climates. Hence it is im- 
portant to understand their form and origin 0]. The 
large number and variety of geological processes acting 
in the formation of landscapes makes this a problem of 
formidable complexity. Numerous hypothetical models 
have been proposed to help understand phenomena com- 
mon to many river systems on Earth P, H, p[. In this 
paper, we approach to problem by studying fluid flow 
on continuous surfaces. This straightforward approach 
leads to insights about the origin of observed statistical 
features in real landscapes. 

Much attention has been devoted to the so-called 
slope-area law, that presumably relates the erosion of a 
landscape to the flow of water on it |p, O, R, H, H, nfl, 1^, 



[l^ , [13 , |14| . The area, projected onto the horizontal, that 
discharges its rain water into a surface element downhill is 
referred to as the drainage area a of the surface element. 
The observed correlation of local slopes with drainage 
area has often been attributed to erosion mechanisms [0 . 
A popular hypothesis is to balance downward sediment 
mass transport with tectonic uplift (e.g. [|l^ |l^, p^). 
Here we show that such a correlation is found on simple 
surfaces for purely geometric reasons. The main results 
of this paper are an exact relation, eq. (^, between slope 
and convergence for a wide class of random landscapes, 
and a geometrical explanation of the slope-area relation. 
The basic, intuitive reason for both results is that flatter 
points on a surface have a stronger tendency to accu- 
mulate water from their neighborhood. Flow aggregates 
preferably on flatter regions. 

The primary distinction between our approach and 
that of many previous investigations is that we look at 
the geometrical effects of a surface rather than the active 
role of rivers. Landscapes can be actively formed by con- 



tinual erosion, but where do rivers choose to form on a 
prescribed landscape? What are their statistical proper- 
ties in the absence of any physical mechanism altering the 
landscape? Studying flow on prescribed surfaces allows 
to understand the passive role of rivers. 

A useful family of artificial landscapes are Gaussian 
surfaces, which are formally introduced in the text. They 
form a broad, yet relatively simple class of random sur- 
faces. Their topography has rich implications upon the 
statistical properties of flow on them. Gaussian surfaces 
may therefore serve as a "null hypothesis" that repre- 
sents basic geometric effects. To determine the traces of 
erosion in landscape topography, one needs to find out 
what is special in the statistics beyond these geometrical 
effects. 

In Section II we ask how the geometric form of a 
surface determines the drainage behavior, leading to a 
differential expression for the topographic convergence 
of a surface. Section III deals with Gaussian surfaces, 
which serve as our model system. In the fourth section 
we demonstrate that a slope- area relation is present on 
Gaussian surfaces and its explanation applies also to real 
landscapes. In Section ^ the slope-area relation is ex- 
plored more quantitatively and its limitations are pointed 
out. Conclusions are reiterated and discussed in the last 
section. 



II. FLOW CONVERGENCE ON A SURFACE 

We consider a smooth surface described by its height 
h{x,y) as a function of the horizontal coordinates x and 
y. Water is assumed to follow the path of steepest de- 
scent, neglecting its inertia and volume. 

To begin, we define and calculate the convergence of 
flow on the surface. Intuitively, a flow is convergent when 
neighboring fiow lines approach each other as they con- 
tinue downhill. Consider the paths of two flow lines sep- 



arated by an infinitesimally small distance t at height h. 
A little further downhill, at a height h — 5, they are sepa- 
rated by a distance t'. Figure || illustrates the situation. 
The rate of lateral approach is given by 



P 



\t\-W\ 



(1) 



meaning that the convergence p is the relative contraction 
of the contour line segment per unit height. Dividing by S 
makes the definition independent of the spacing between 
contours. 

To obtain p, we first find the direction of the tangent to 
the contour, described hyt = {x{t), y{t)). The parameter 
t can be chosen to represent the length of the contour 
element. We denote h{t) = h{x{t),y{t)). A first-order 
expansion 

h{x, y) = /i(0, 0) + 9:^/1(0, 0)x + dyh{Q, 0)y + 0{x^ + y^) 

(2) 
yields the equation for the tangent t of the contour at 
point (0,0). Since h{t) — const, and It] = t, one obtains 
to first order 



dyh d^h 
t , , and y = —t- -. 

\\Jh\ ^ |V/i| 



(3) 



Derivatives without their argument are here always un- 
derstood to be taken at the origin, that is, we write dxh 
in lieu of d^hi^^^"), and so on. The first-order expansion 
of the gradient along the contour becomes 



V/i(i) = V/i(0)' 



t 



dxxhdy h - dxy hdx h 

|V/i(0)| V dxyhdyh - dyyhdxh 



-0{f 



According to Fig. |l|, 

,7 _ r , ^ V/i(t) 



(4) 



t' = t 



V/i(0) 



\S/h{t)\ |V/i(t)| |V/i(0)| |V/i(0)| 

= i*+ ^T^^yp [v/i(t) - v/i(0)] + o{t5). (5) 

A straight-forward calculation using eqs. (|l|), (Q), and (||) 
then yields 

_ {dyhfdxxh - 2(dxh){dyh)dxyh + {dxh)'^dyyh 
^ " iidxh)^ + idyh)^]^ 

The surface is convergent when the expression (0) is pos- 
itive and divergent when it is negative. The expression 
for the curvature of the contour differs from eq. (H) only 
by a denominator of |V/ip jl^l instead of |V/i|^, so that 
the convergence may also be written as the curvature of 
the contour, k, divided by the slope, p = k/|V/i|. 

The definition of convergence (El) is merely based on 
the geometry of the flow paths and does not involve ve- 
locities. In fact, without further assumptions one has no 
knowledge of velocities and fluxes on the surface. 

Essentially the same notion of convergence as in eq. d3) 
has been proposed earlier |17| and employed in stud- 
ies of landscapes [f^ 11^. In rcf. |17 landscapes were 



decomposed into elements by tracing their height con- 
tours. The generated contour segments were then used 
to evaluate (|il — |i'|)/(|il + \t'\) for a given 6. For an in- 
finitesimally fine grid this agrees, up to a prefactor, with 
eqs. (P and (^. Others (see refs. |g, |l^) have proposed 
the Laplacian W^h as a local definition of convergence, 
which can differ from eq. (m even in its sign. The use 
of the Laplacian applies for a particular model of ffuxes 
and erosion, but it is not a general expression for to- 
pographic convergence. For example, the hyperbolic sur- 
face h{x, y) = 2x^ — y^ would be erroneously identified as 
convergent everywhere. The proper distinction between 
topographic convergence and divergence is given by the 
sign of eq. (||). 

It is worthwhile to pause here to contemplate the intu- 
itive reasons for convergence. If the contours are strongly 
curved the surface acts like a funnel and water converges 
rapidly. The slope influences the convergence via the 
curvature, but it also alters the available path length per 
unit height. The latter effect is reflected in the denomi- 
nator oip = n/\Vh\. Fig. demonstrates this slope effect 
in a simple example. In this example, both k and |V/i| 
change to increase the convergence. On steep slopes wa- 
ter falls quickly along paths with little relative motion, 
while on ffatter slopes paths are deffected more easily. 

The notion of convergence intrinsically requires that 
the elevation is a function of two spatial coordinates. It 
has no equivalent in a one-dimensional setting, where the 
height is considered as a function of one coordinate only. 
Hence, any consequences derived from the convergence 
expression are three-dimensional in nature (height and 
two horizontal coordinates). 

In our idealized setting rivers may be associated with 
locations carrying much more running water than neigh- 
boring points at the same height. They can do this by 
accumulating water from uphill more effectively. If a 
point converges flow faster than a neighboring point at 
the same height, fluid moves towards the point of higher 
convergence. Consequently, locations of strongest con- 
vergence tend to act as rivers. It is thus clear that rivers 
tend to form in locations of large p, where the slope is 
low and the curvature high. Rivers are thus associated 
with a local optimality property, resulting from the mere 
fact that they flow on surfaces. This optimality property 
is not directly equivalent to "optimal channel networks" 
[Q or any of the other proposed optimality principles |2^ 
for river networks. 



III. GAUSSIAN SURFACES 

Gaussian surfaces are random continuous surfaces de- 
scribed by a superposition of waves with random phases 
pl| , p^ , |2^ , p4[ . Results on Gaussian surfaces have been 
applied in various fields, including electric noise |2^, wa- 
ter waves 1^, |3|, beaches |2^, image processing [ p7[ , 
and the clustering of galaxies |gg. We will here need 
only some of their elementary properties. 



A Gaussian surface is defined as 

h{x, v)^^ c{k, ly-^ik^+^y+vi'^^'!-)) 



(7) 



where (^ is a random variable uniformly distributed be- 
tween and l-K. The amplitudes c{k,£) are arbitrary, 
but, of course, they need to decay sufficiently fast with 
wavenumber to guarantee the continuity and differentia- 
bility of the surface. Any Gaussian surface has the prop- 
erty that its heights and all of its derivatives are dis- 
tributed Gaussian. Furthermore, the first derivative is 
statistically independent of the height, and the second 
derivatives are independent of the first derivatives |21| . 

This last property allows us to evaluate the conver- 
gence for a given slope quantitatively. Qualitatively, it is 
already clear from our early discussion of expression IM) 
that flatter areas tend to converge flow more rapidly. For 
Gaussian surfaces this relation can be made precise. The 
convergence, p, averaged over all points of a given slope, 
|V/i|, is called the conditional average of the convergence 
on the slope, (p| |V/i|). The simple conditional aver- 
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age (p|(V/i)^) = vanishes, because p changes sign as 
h -^ ~h. Thus we look at the square of the convergence: 



+ {{dyyhf) 



^A{{d.,yhf) 



idyhY 






{Vhf 



(vhY 



my 



{d^h)\dyhf 



(V/i)8 



{Vhf 



+2 {{dxxh)dyyh) 



{dxh)^{dyhf 



{vhr 



{vhy 



Here, we have made use of the statistical independence 
of first and second derivatives. The mixed terms of p^ 
vanish because of the statistical symmetries x — > —x 
and y ^ —y. The conditional averages now reduce to 
geometrical factors, which can be rewritten for a sta- 
tistically isotropic situation in terms of an angle 9 = 
ta,n^^ (dyh/dxh): 



(p^Vhr) = {id^^^hr){sin^ 6) ^^+{idyyhf ){cos^e)^^ 



+4:{{dxyh)^){sin^dcos^e) 



1 



(vhy 



-2 {{dxxh)dyyh) {sm^Ocos^ 9) 



(Vhy 



I {{dxxhf) + \ {{dxyhf) + \ {{dxxh)dyyh) + ^ {{dyyhf) 

l/{Vhf 



(vhy 



(8) 

(9) 
(10) 



r 



Therefore, the convergence behaves as the inverse of slope 
squared. 



\p\ - l/i^hf 



(11) 



or vice versa, |V/i[ ~ 1/a/|p|. Since the slope enters 
quadratically in ( [ll| ) , it has a strong efii'ect on the conver- 
gence. This is the relation between convergence and slope 
for a Gaussian surface, irrespective of the surface's spatial 
correlations. In fact, for result (pi) to be valid, the sur- 
face does not necessarily need to be Gaussian. Instead, 
one requires only that the first and second derivatives of 
height be statistically independent of each other. 

We can check this result on both synthetic and real 
landscapes. A Gaussian surface can be efficiently cre- 
ated using Fourier transforms |29 . Interpolation could 



be used to obtain heights and gradients at any position 
between grid points. Real landscapes from digital eleva- 
tion maps, which we use below, are however convention- 
ally evaluated on a discrete grid. We opt to evaluate syn- 



thetic surfaces and real landscapes in the same discrete 
way. Fluid paths are determined by following the direc- 
tion of steepest descent among the eight neighbors. Fluid 
trajectories are stopped when they reach a local mini- 
mum, implying a lake or ocean at that place. We choose 
the smallest wavelength A as four grid cells long, which is 
twice the smallest representable wavelength. The largest 
wavelength A equals the domain size. 

Figure shows the convergence-slope relation mea- 
sured in the simulated surfaces, which must agree with 
formula (pi). The filled dots are the measurements from 
a simulated surface and the dashed line is the theoretical 
calculation using the prefactor in eq. (0). The tiny de- 
viations between theory and numerics have two sources: 
The average was only taken over one realization and the 
domain is periodic in both directions. In the same fig- 
ure are data from real landscapes, discussed later in the 
text. Note now, however, the good fit over two orders of 
magnitude in slope; almost the entire range. 



The result (^ shown in Fig. || is independent of the 
choice for the power spectrum. Yet, for completeness, 
we note that the spectrum used in this example is a 
power-law spectrum, where the coefficients in eq. ([7|) are 
\c{k)\ oc l/|fcp+''^ and the phases are random. The 
exponent p is known as the Hurst exponent. It has 
a simple relation with the "roughness exponent" that 
characterizes the spatial correlation of the height. The 
structure function (\h{x) ^ h{x + R)\'^') is proportional 
to R^P for < p < 1. If p is small, one has to intro- 
duce a cut-off scale to maintain the differentiability of 
the surface, and therefore we take \c{k)\ = l/|fcp+''' for 
2tt/\ < k < 2tt/A and c{k) — otherwise. Graphs with 
solid dots, in Fig. raand elsewhere, correspond to a Hurst 
exponent of p = 0.7. 

In concluding this section, we note that convergence 
is intuitively related to the drainage area. Points of 
stronger convergence will, on average, possess larger 
drainage areas, because they collect more flow than 
points of low convergence. If p increases so does a, lead- 
ing to a slope-area relation. 



IV. SLOPE-AREA RELATION 

Erosion rates increase with the amount of flowing water 
and are also affected by the local slope of the landscape, 
among other factors |2|. If other parameters, like width, 
depth, and velocity of a river are considered a function 
of total discharge and slope, one arrives at the conven- 
tional assumption that erosion is a function of slope and 
drainage area only H. Local topography can carry the 
signature of the entire basin draining into it ||30|] . 

On real landscapes a statistical relation is observed be- 
tween the local slope and the drainage area a. Its quan- 
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titative form is roughly \Vh\ ^ 1/V^ 
are however substantial variations of the exponent of 0.5 
for different river systems [pi. We have analyzed real 
landscapes from digital elevation maps of coastal land- 
scapes from the Mendocino region in northern California 
| p2| , 33 . A detailed description of the Juan Creek data 
is given in ref. p2| and the Noyo river basin is discussed 
in ref. 1 33 . These data have 30m resolution in both hor- 
izontal directions. Fig. ^a shows the slope-area relation 
of a Gaussian surface together with that of three real 
landscapes. One sees that the random surface exhibits 
a slope-area relation similar to the one seen in the real 



data l34 



For Fig. ]m, the slope was calculated as the maximum 
drop from pixel to pixel. Another method, used in all 
other flgures, is to take a finite difference in both hor- 
izontal directions and to evaluate the absolute value of 
the gradient in this way. Using a stencil of five points 
in both direction, this is two orders more accurate than 
the previous method, but not necessarily more realistic, 
because the gradient is in this case not exactly aligned 
with the direction of the flow. Figure Hb repeats Fig. Ha, 



with the alternative slope measurement. The observed 
slope-area relationships are weaker in this case, but still 
roughly equivalent for real and Gaussian surfaces. The 
discrepancy with Fig. ^a arises merely from the ambigu- 
ity of measuring the slopes from the data [ p5| . 

Real landscapes are more complex than Gaussian sur- 
faces. Erosion might alter landscapes in a way not cap- 
tured by Gaussian surfaces. Yet, the direct compar- 
isons of the slope-area relation in Fig. ^ shows a sim- 
ilarity between them. The extent to which the geo- 
metric explanation of the slope-area relation applies to 
real landscapes can be further probed by looking at the 
convergence-slope and convergence- area relation of the 
real landscapes. Fig. ^ shows that the convergence-slope 
relation of Gaussian surfaces is reasonably true also for 
real landscapes. All of them show a behavior close to 
p^ ^ l/lV/ij"*. Fig. shows the convergence-area rela- 
tion for real landscapes together with that for a synthetic 
landscape. Particularly small drainage areas are found 
for divergent points. The minimum drainage is the size 
of one cell, which is 900m^. After a rapid initial growth 
of the average convergence for smaller areas, the graph 
approaches its asymptotic behavior. Again, we see that 
real landscapes behave in this respect similar to Gaussian 
surfaces. Hence, there is no doubt that the geometric ef- 
fect found for Gaussian landscapes is also present in real 
landscapes. 



V. VARIATIONS IN THE SLOPE-AREA 
RELATION ON GAUSSIAN SURFACES 

The slope-area relation depends on the statistical pa- 
rameters of the surface, unlike eq. (J^). Therefore, we ob- 
tain pertinent statistics from numerical simulations. In 
the following, we discuss how the relation changes with 
the choice of the power spectrum. We explore various 
surfaces that vary over a wide range of length-scales. 

Fig. g shows the slope-area relation for power laws 
with several different Hurst exponents p, one exponential 
power spectrum, and a landscape that is half-flUed with 
water (ignoring the lower half when taking the statistics). 
For the surfaces with short spatial correlation there is a 
clear area dependence of the slope, while the surfaces that 
change more slowly exhibit less of a dependence. For 
landscapes that are statistically self-similar, and hence 
characterized by their Hurst exponent p, the slope-area 
relationship becomes more and more pronounced as p 
goes from one to zero. 

Real landscapes are said to have Hurst exponents of 
about 0.7 ^2l: although there are considerable variations 
among individual landscapes (see p% and refs. in p8[). 
We determined the structure functions of the real land- 
scapes used in this study. They exhibit only a rather 
limited range of scaling, if at all, but they tend to be sub- 
stantially flatter than 0.7. This corresponds to a lower 
Hurst exponent and thus a faster decorrelation. However, 
given the average tilting of the coastal basins and their 



small size, this might not be significant. In any case, in 
the range of Hurst exponents < p < 1 the simulated 
Gaussian surfaces show a clear slope-area dependence, as 
seen in Fig. ^. The graphs for p — OA and p — 0.7 have a 
slope in the range of 0.15-0.25. Varying the power spec- 
trum can lead to flatter or steeper slopes, but we have 
never observed a slope-area graph steeper than about 0.3. 

Fig. ^ shows that for large convergence, p grows 
slower than the area a for the real and the simulated 
landscapes. The same is true for simulated landscapes 
with other parameters (not shown in Fig. H) and we have 
not found any Gaussian landscape that would violate 
this trend. Since p < 0{a) and p ~ l/|Vft.p, it follows 
that |V/i| > 0{l/y/a). Hence, one expects that slope 
decays no faster than \j ^fa. As mentioned above, the 
actual restriction appears to be even stronger, around 
^-0.3^ This indicates the presence of a bound for the 
area-dependence of the slope on a Gaussian landscape. 
The convergence depends strongly on the slope and the 
area depends strongly on the convergence, and therefore 
the area strongly on the slope, or, equivalently, the slope 
weakly on the area. Hence, it is difficult to get a par- 
ticularly strong slope-area relation through geometrical 
effects alone. 

In the particular case of Fig. ^ the graphs for the real 
landscapes change their slope at about lO^m^. This cor- 
responds geologically to a transition from coUuvial (loose 
debris) channels to bedrock channels [^. The fact that 
the simulations nowhere exhibit the final drop in the 
slope-area relation to its full extent seen in real land- 
scapes (Fig. H), indicates that there are contributions 
to the slope-area relation beyond the geometrical effect 
we have unearthed in this paper. Exceptional cases of 
real landscapes are known with slope-area relations that 
are substantially more pronounced than the ones we en- 
counter on Gaussian surfaces [Q ^. This leaves room 
for erosion, tectonics, or further mechanisms contributing 
to the slope-area relation. However, since a substantial 
slope-area relation is produced by Gaussian surfaces with 
realistic Hurst exponents, the slope-area relation carries 
less information about these processes than previously 
thought. 



VI. DISCUSSION 

In this work, we found and explained a correlation be- 
tween slope and drainage-basin area on random topogra- 
phy. We gave a two-step argument. The first step was 
to derive analytically the convergence expression and to 
show explicitly the dependence of convergence on slope. 
The second step connected the convergence with the area 
of drainage. This connection is intuitive, but we also 
demonstrated it numerically. Further, we have left lit- 
tle doubt that the same effect appears on real landscapes 
too, by verifying the validity of both steps with real data. 



We have also found that the slope-area correlation 
cannot be strengthened indefinitely by varying the pa- 
rameters. The slope-area relation is often approximated 
as |V/i| ^ a~^ . Since we find for random surfaces 
< < 0.3, an exponent in this range implies nothing 
about geological processes. It is instead only a conse- 
quence of water fiowing downhill on random landscapes. 
This is in direct contrast to hypotheses that attribute any 
slope-area relation to erosional mechanisms (see refs. in 

i). 

If erosional mechanisms were responsible for slope- area 
relations one would expect a deterministic relation be- 
tween slope and area. A certain drainage area would 
always lead to the same slope at its outlet. Fluctuations 
would then be due to variations in geological parameters. 
On Gaussian surfaces on the other hand there are intrin- 
sic statistical fluctuations and the slope-area relation is 
a statistical correlation rather than a deterministic for- 
mula. Investigations of qualitative differences between 
slope-area fluctuations on real and random landscapes 
are therefore expected to yield additional insight. 

Although this study is motivated by geomorphology 
the same statistical findings are applicable to other phys- 
ical systems, where the motion is strictly along the gra- 
dient of the force field and the inertia negligible. The 
slope-area relation in this context states that for a quickly 
decorrelating force field the domain of attraction will, on 
average, be larger when the force is smaller. (Ref. |24] re- 
views applications of Gaussian surfaces when the motion 
is along equi-potential lines, rather than perpendicular 
to them.) 

The slope-are relation is but one consequence of the 
convergence expression y^. Another consequence, as we 
have briefly discussed, is the formation of rivers prefer- 
ably on locations of high convergence. The local opti- 
mum is with respect to points on the same contour line, 
that is, p is maximum along h = const., or in terms of 
our earlier notation dp/dt = and d^p/dt^ < 0. This 
indicates a previously overlooked optimality principle for 
the organization of rivers. (Note that this extremal prin- 
ciple is not strictly valid for points at different heights. 
If it were, it would immediately imply that the longi- 
tudinal stream profiles are concave). Whether such an 
observation may by of any use is unclear. What is clear, 
however, is that Gaussian or otherwise random topog- 
raphy has far-reaching non-trivial implications on basin 
statistics. 
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FIG. 1: Schematic drawing for the calculation of flow conver- 
gence. The dotted curved hnes are height contours. 
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FIG. 2: Illustration of the effect of slope on the convergence 
of flow lines. A simple surface of the form h{x, y) — x^ + by 
is used as an easily visualized example. The solid lines in the 
upper panels of the figure show the flow lines for (a) 6 = 1 and 
(b) & = 3. The dotted lines are height contours. The contour 
interval is the same in both parts of the figure. One sees the 
stronger convergence in the flatter of the two surfaces, (a) , as 
well as the stronger convergence in the flatter regions towards 
the center of the x-axis within each flgure. The behavior of 
the convergence p{x) and the slope s{x) (actually s{x)/2) are 
illustrated in the lower panels. 
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FIG. 3: Convergence-slope relation for synthetic (filled dots) 
and real landscapes (unfilled symbols) . The dashed line shows 
the analytical result for the synthetic landscape. All data 
are averaged over logarithmically spaced bins on the abscissa. 
The unfilled symbols correspond to three different landscapes 
in Northern California as described in the text. (The units 
for the synthetic landscape are arbitrary and hence only the 
slope of the graph can be compared to real data. The two 
lowest graphs are shifted downwards by factors of 40 and 1600, 
respectively, for better visibility.) 
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FIG. 4: Slope-area relation for synthetic and real landscapes. 
The solid line with filled dots in parts (a) and (b) is for a sim- 
ulated Gaussian surface at a resolution of 4096 x 4096. The 
length units for the simulated landscape are chosen arbitrarily. 
The unfilled symbols correspond to Juan Creek {<>), Middle 
fork of the Noyo river (o) , and the North fork of the Noyo river 
(A). These symbols are used uniformly in all figures. (The 
graphs for the last two landscapes are shifted downwards by 
factors of 1.4 and 3, respectively, for better visibility.) The 
two parts of the figure differ only in the way the slope is de- 
termined. In both parts real and artificial landscapes exhibit 
similar slope-area graphs, suggesting the relation has a simple 
geometric origin. 
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FIG. 5; Convergence-area relation for synthetic (filled sym- 
bols) and real landscapes (empty symbols). The two plots 
show the same data, (p|a), on (a) a semi-logarithmic and (b) 
a log-log plot. The length units for the simulated landscape 
are chosen arbitrarily. The dashed line in (b) corresponds to 
p (X a. 



Q. 
O 

CO 



10' 



10" 



10" 




10" 



-q— 1 r-q— 1 |-q— i i-q- 



' - exp. 
p=0.4 
p=0.7 w/ ocean 

p=0.7 
P=1 

p=1.4 . 




jj I I—lI I I— J I I— lI I L-J I I— lL 



10^ 10"^ 10^ 10^ 10'' 10^ 10^ 



FIG. 6: Slope-area relation for Gaussian landscapes. The dif- 
ferent graphs correspond to landscapes with power-law spec- 
tra with different Hurst exponents, an exponential power 
spectrum, and a landscape half-filled with an "ocean" . They 
are obtained at a resolution of 1024 x 1024. The length units 
are chosen arbitrarily. 



